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I. INTRODUCTION 


In the search for an almost infinite and clean energy source mankind has looked 
towards nuclear fusion. In this endeavor the understanding of the properties of plasmas 
is a necessarv undertaking. Since we are unable to produce machines of infinite size, 
we must pay particular attention to the interactions between the plasma and the surfaces 
of the machines used to contain it. 

The plasma-surface interaction occurs primarily through the following erosion 
processes: 

e Sputtering 

e Unipolar Arcing 

e Heat pulses 

¢ Gases which implant into the surface can generate chemical erosion, embrittlement 

and blistering. 

These processes are of concern due to their limiting the life of equipment components 
and their resultant introduction of impurities into the plasma. Of these processes it 1s 
widely believed that Unipolar Arcing 1s dominant. 

Most research into Unipolar Arcing has thus far investigated onlv the plasma above 
the surface. This thesis starts an investigation of the dynamics of the surface itself. It 
begins with experiments to determine if the bulk emission of material has any directional 
dependance. The results of these experiments indicate that not onlv is there directional 
dependance but there is also more structure to the ejected material than previously 
thought. 

Thin films are used for two purposes: to see how the same arc interacts with two 
different materials and to find if there is any correlation between pit depth and diameter. 
It was found that the amount of damage may be dependant upon the material’s melting 
and boiling points. This also lead to a method for calculating the arc current and to the 
role of surface tension. Investigation into the role of the surface tension lead to the re- 
alization of the similarity between Unipolar Arcing and the physics of arc welding. This 
relationship 1s exploited in Chapter V which provides a preliminary theory for the sur- 
face dynamics. No correlation was found between pit depth and diameter, however 


these experiments indicate that, for the parameters used, the maximum pit depth is be- 


tween 5.5 and 10 microns. In all this research the Unipolar Arcs are generated by a 


plasma which itself is produced by a CO2 laser. 


td 


Il. BACKGROUND AND THEORY 


How high energy laser pulses and plasmas interact with various material surfaces 
has long been under study. Of particular interest has been how laser pulses damage a 
target surface when the energy of the photons in those pulses is less than the energy 
necessary for photo-ionization of the material. Presented here is one of the current 


theories describing this process, it involves the process known as unipolar arcing. 


A. BACKGROUND 

Before going into the theories delving into unipolar arcing some background must 
be given on plasmas and their interactions with surfaces. This section begins with a 
definition of what a plasma is and mentions some of its properties. The Debve length 
is derived and a brief description 1s given on how a plasma 1s produced by a pulse from 
a CO2 laser. 

1. Definition of a plasina 

In nature there are four states of matter, these are solid, liquid. gas and plasma. 

These states represent increasing levels of disorder and require higher energy levels with 
progress from one to the next. The solid is the state of lowest energy and consists of 2 
concentration of atoms or molecules bound together in a rigid lattice structure by their 
mutual electrical attraction. As energy, heat, is added these atoms will vibrate more and 
more strongly until they can no longer maintain their lattice structure. In this regime 
the atoms can move about but their energy is not sufficient to fully overcome their mu- 
tual attraction. This is the liquid phase. When more energy is added so that the atoms 
can completely overcome their mutual electrical attraction they enter the gaseous state. 
To form a plasma it 1s necessary to ionize a portion of this gas. That is, sufficient energy 
must be given to at least one electron in an atom to free it from that atom. This results 
in a gas composed of both neutral and charged particles. Chen [Ref. 1, p.3] defines a 


plasma as 


A plasma 1s a quasineutral gas of charged and neutral particles which exhibit col- 
lective behavior. 


The term ‘quasineutral’ simply means that the gas contains approximately equal con- 


centrations of ions, 7, and electrons, ”,. such that 


Dd ed {1} 


where n is called the plasma density. Collective behavior comes from the charged par- 
ticles. These particles produce and are affected by electric and magnetic fields. Because 
of this, elements of the plasma will exert forces on one another even at a considerable 
distance. Thus they will act together, hence the term collective behavior. 

Temperature is another property of a plasma. If it is in thermal equilibrium 
then the temperature is related to the average kinetic energy of its particles by the 


E joo 2 


Here k 1s the Boltzmann constant 
k= 1.38 x 107? WK {3} 


and T is the absolute temperature in kelvins. 
2. Debye shielding 

Debve shielding refers to the ability of a plasma to shield out electric potentials 
which are apphed to it. When an electric field 1s applied to a plasma the electrons and 
ions will drift in opposite directions in response to this field. The charge separation be- 
tween the electrons and ions acts to counter the applied electric field until it 1s essentially 
cancelled out. For instance if a positively charged ball is placed in the plasma it will 
attract the electrons and repel the ions. The electrons will form a cloud of equal charge 
around the ball thus effectively cancelling its electric field. 

An approximate width of the charge separation caused bv a potential, , can 
be easily deduced if all of the motion 1s considered to be done by the electrons. The fact 
that the 1ons are so much more ma.sive than the electrons makes this a good assur p- 
tion. The ions will then form a uniform background upon which the electron motion 


takes place. The ion density will be the same locally as it 1s at a large distance, hence 
HM; — co {4} 


Assuming that the potential is held at a fixed value, ®,, on a plane at x=0, then D(x) 


is desired. Poisson’s equation in one dimension is 


2 
PEG ee c(i). (5} 
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The Boltzmann relation for the density of electrons 1s 








e® 
Ne = Nee kT, {6} 
so that equation (5) becomes 
2 eD 
E,- L = eNge| € kT, —I1 |]. {7} 
ax 
The solution to this equation 1s 
aes 
D(x) =D,e 4p {$} 


where /, is the debve length and 1s defined as 


er lele. 
i) a mame {9} 
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The debve length 1s a measure of the shielding distance. It is the distance over which the 
potential is reduced to l’e of its maximum value. 
3. Plasma generation by the CO2 laser 
In this research the plasina was generated by a CO2 laser. The laser sends an 
electromagnetic pulse to the target. If the target is considered to be at x=0 and the 
pulse travelling to the left on the x-axis, then the electric field of the incoming laser 


meise, E,, is given by 





oS sin( ee tor) {10} 


maximum amplitude of the E-field 


where E,, 


A= laser wavelength 
w= angular frequency 
[= (indie. 


For simplicity the target is assumed to be perfectly reflecting. The electric field of the 


reflected pulse is then 





Ev=-E a sin( ns -wtte] {1 


where € 1s a phase angle. Perfect reflection also provides that 
ele ee leo. {13} 


If the target 1s a perfect conductor, then the electric field at x=0 parallel to its surface 
must be zero. Thus ¢=0. The superposition of these waves produces a standing wave 


whose electric field is given bv 





ee 2 ax 
E=2E, sin( 5 ) cos(wt) {14} 
and whose intensitv 1s 
2 22 ees 
[Ie = 45 sin ( 255). {15} 


Assuming that an inital electron concentration exists, then electrons above the surface 
Will oscillate in response to this electric field. Maximum oscillations will occur at the 
maxima of the standing wave, 1.e. at x = 4/4, 34/4, etc... The source of this initial 
electron concentrauon Is as vet unclear. 

An actual target surface is not perfectly reflecting, thus some of the energy of 
the laser pulse will be absorbed bv that surface. As the target is heated neutral particles 
adsorbed will be released. Oscillating electrons will collide with these neutrals ionizing 
manv of them. When thev are ionized the neutrals release electrons which in turn os- 
cillate and collide producing further ionizations. This process continues until the plasma 
density rises to a point where the plasma frequency is equal to the laser frequency. Thie 


laser will then be reflected by the plasma. The plasma frequency, w,, 1s given by 





es |= {16} 


where e and m are the charge and mass of the electron. Thus the plasma density ob- 


tained will be 


2 
EIN 


Aa=—— {17} 
me 


which for the CO2 laser (/ = 10.6 micron) vields a critical plasma density of 
n=9.93 x 10% m7, 


B. THE ROBSON-THONEMAN MODEL 

In 1958 Robson and Thoneman [Ref. 2] observed what they described to be unipolar 
arcing occurring on a mercury electrode in the presence of a high energy plasma. Their 
explanation of this phenomenon follows. 

Whenever a metal plate 1s exposed to a plasma it will collect both the positive ions 
and the negative electrons. The electrons having much higher random velocities than 
the ions will initially collect on the plate more rapidly. An electrically isolated metal 
plate will collect a negative charge and thus a negative potential. This negative potential 
will repel electrons even as the now more positive potential of the plasma near the sur- 
face attracts them. Onlv the electrons at the high end of their maxwellian distribution 
will have sufficient energy to overcome these potentials and reach the plate. At equi- 
librium this potential is known as the floating potential (Figure 1). 


The negatively charged plate receives an ion current density of 


Ji = Nie; {19} 


Where »;=10n density 
é; =10n charge 
;= mean velocity with which the ions reach the sheath boundry 


y 
silos 

a ( M 
Af = ion mass 


7, = electron temperature 
o = factor on the order of unity. 





Given that the plate potential with respect to the plasma is V, then the electron current 
density is 


eV 
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Figure 1. Equilibrium Plasma-Surface Interaction Prior to Cathode Spot Forma- 


tion. 


where 7, = electron density 
e, = electron charge 
¥. = mean random Velocity of.the plasma electrons 


ri SkT, : 
= THN 


m= electron mass. 


In steady state the plate is at the floating potential, 1’, and the total current 1s 








Jit ea {23} 
Assuming that a1 and e, = — e, (1.e. singly charged ions) then we have 
Jt ar) 
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This equation shows that given a sufficiently high electron temperature, 7,, the floating 
potential, V,, will surpass that needed to sustain an arc. Such an arc current of electrons 
from the plate to the plasma will lower the local potential between the plate and plasma 
to the cathode fall potential of the arc, V,. This lower potential will allow more electrons 
from the high energy end of their maxwellian distribution to cross the sheath from the 
plasma to the plate thus completing the circuit (Figure 2). There is then a circulating 
current which maintains the requirement that the total current to the plate be zero. The 


circulating current is described by 


KY, E,U eV 
wii e kT, —ekfT, 
£ ‘ 2mm 4 


Where A = plate area exposed to the plasma. Thus this is indeed a unipolar arc, an arc 





which requires only one electrode and runs off the thermal energy of the plasma 
electrons. 

Robson and Thoneman tested their theory by creating a plasma in a vacuum tube 
and then heating it with RF energy. In their experiments they found good evidence in 


support of their theory and observed unipolar arcing taking place. 
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Figure 2. Plasma-Surface Interaction During Unipolar Arcing. 


C. THE SCHWIRZKE MODEL 

In 1980 Schwirzke [Ref. 3] expanded upon the Robson-Thoneman model, he con- 
nected the unipolar arc phenomenon with the laser-target interaction. A Q-switched 
neodymium laser was used in his experiments. Later research by Schwirzke found gen- 
eral applicability to Tokamaks [Ref. 4] and the CO2 TEA laser. A plasma in contact 


with a surface will form a sheath with the floating potential given by 


eee i 
v= ye in( Sey ) {26} 
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It was also stipulated that the sheath width is proportional to the Debye length 


fee 
An=, | oe {27} 
Snne~ 


Thus the electric field across this sheath and 1n a direction normal to the surface will be 
given by 


, 


} 25 eee 
\E |~——=.fankT, n( us ) {28} 
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Umipolar arcing develops if the sheath potential is high enough to sustain an arc. The 
electric field is then intense enough to cause the emission of electrons from a surface spot 
into the plasma. “The minimum requirement for the onset of arcing is that the sheath 
potential is comparable to the ionization energy ...” [Ref. 3] of the material. 

“Por an arc to develop it is also necessary that the ion density increases in front of 
the cathode spot in order that a larger electron current can flow into the plasma.”[Ref. 
3} This occurs when the electric field concentrates around protrusions or ‘whiskers’ 
(Figure 3) causing increased ion bombardment and hence heating. The increased local 
heating causes increased desorption of surface gases, oils and metal evaporation. This 
causes an increased concentr~tion of neutrals in the sheath of which a small percentage 
is ionized. The small percentage of ionized neutrals causes a drastic increase in the local 
ion concentration. These new ions are then accelerated by the sheath potential back to 
the surface causing further increases in surface temperature, vaporization and electron 
emission from the hot cathode spot is the result. The locally increased 10n concentration 


above the cathode spot leads to a radial electric field (Figure 4) given by 
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Figure 3. Local Electric Field Enhancement at a Whisker Tip. 
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“This reduces the potential in a ring-like area surrounding the higher pressure region 





above the arc spot. The sheath potential will be lowered there and more electrons from 
the high energy tail of the Maxwellian distribution can reach the surface, thus closing 


the current loop.” [Ref. 3] 


D. THIN FILMS 
Thin films [Ref. 5] have been prepared for quite a long time. Thin film technology 

generally refers to the deposition of layers by an atom by atom process as opposed to 
lavers formed by laying down relatively large particles. One of the most inportant aims 
in making a thin film is generating continuous layers, large particles may or may not 
react together to form such layers. The basic techniques for the deposition of thin films 
are: 

e Electroplating 

e Chemical Reduction 

¢ Vacuum Evaporation 


°¢ Sputtering 
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Figure 4. Schwirzke Unipolar Arc Model. 


¢ Chemical Vapor Deposition = _ 


In the research reported in this thesis the vacuum evaporation method was used, The 
method of thin film production by evaporation consists of four physical stages: 

¢ conversion of the material to be deposited into a gaseous state; 

® transport to the substrate; 

e deposition onto the substrate; 


® rearrangement or modifications of their binding as more material arrives.[Ref. 6] 


Vac’ um evaporation requires a relatively clean vacuum in a range typically obtained 
by a modern ion or diffusion pump system. In this vacuum a vapor source will develop 


a rate of vaporization given by 


Se 
Neos oe 1022 Qiie molecules (30} 


JM,T — cm*sec 


Where a, = evaporation coefficient (+1 for a clean evaporant source) 
AM, = element atomic weight of the vapor source 
P, = vapor pressure of the vapor source 
7 = temperature. 


How this vapor finally deposits onto the desired substrate is dependent upon the geom- 
etry of the system, that is, it depends upon the geometry of the vapor source, the 
substrate, and relative positions of the two. For example, using the oval dimple boat in 
Figure 5 will give deposition profiles given in Figure 6, whereas using the coil in Figure 


7 will give the deposition profiles shown in Figure 8 [Ref. 7]. 





Figure 5. Oval Dimple Boat. [Ref. 7, p.27] 
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Figure 6. Oval Dimple Boat Evaporation Distribution Profiles. [Ref. 7, p.27] 
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Figure 7. 3/8 inch I.D. X 6 Coil. (Ref. 7, p.110} 
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Figure 8. Coil Evaporation Distribution Profiles. [Ref. 7, p.110] 


E. THIN FILMS AND THE UNIPOLAR ARC 

The purpose of this thesis is to investigate the dynamics of damage on the surface 
caused by unipolar arcing. By gaining an understanding of these dynamics methods of 
preventing damage may be sought after. | 

This being the first investigation dealing exclusively with the dynamics of the sur- 
face, as opposed to the plasma over the surface, the experiments begin simply. While 
the shape of the pit in the surface is as yet unknown, they appear to tunnel normal to 
the surface. If the majority of them are indeed normal to the surface, then the gases they 
eject will be ejected normally. This will occur regardless of the incident angle of the CO2 
laser pulse. The first set of experiments provide a preliminary investigation of this. 

How are different materials affected by unipolar arcing? One method ts to dlumi- 
nate different materials with the laser. While this will provide much data on the gencral 
differences in their interactions, a more interesting experiment is how the same arc af- 
fects two different materials. This the basis of the second set of expermments. In these 
experiments a thin film of one material, copper, is placed over another, tvpe 304 stainless 
steel. This actually serves two purposes. The first is to observe differences in damage 
between the two materials. While it is expected that the film will suffer more damage, 


just how much more damage is the question. The second purpose is to find an estimate 


of the amount of heat generated during the arcing process. Where damage Is restricted 
to the film only, then knowing the material damaged and its volume allows calculation 
of the heat input. Assuming this heat 1s provided by the arc current allows its calcu- 
lation. 

Several assumptions have been made in these experiments. Most have been inter- 
spersed throughout this report and will be brought up as they are needed. The two most 
fundamental assumptions are: 

e¢ the thin films made are of uniform thickness; 


¢ the boundry between the different materials does not significantly affect the arcing. 


Another of the unknowns involved tn the study of the unipolar arc phenomenon 1s 
the depth of the unipolar arc pits. Knowing the depth of these pits will allow measure- 
ment of the actual amount of material ejected from the pit into the surrounding vacuum 
and mav shed further light onto the internal dynamics of the process of the pit forma- 
tion. The final experrment attempts to determine a relationship between the diameter 
of the unrpolar arc pits and their depth. This 1s accomplished using thin films of varving 
thicknesses on a surface which 1s hit by a high power pulse from a CO2 laser. These 
targets are then analvzed to determine the smallest diameter pits which penetrate 
througl. to the substrate for each film thickness. This data 1s then to be used to deter- 


mine a relationship between the size (diameter) of the pits and their depth. 
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Wi. EQUIPMENT AND PROCEDURES 


This research was designed to study the surface dynamics of the molten metal 
produced by unipolar arcing and to find if there 1s a relationship between the diameter 
of the unipolar arc pits and their depth. The equipment used and procedures followed 


to accomplish this objective are described in the following sections. 


A. EQUIPMENT 
The equipment used in the conduct of the research described in this thesis can be 
broken down into three categories: the Vapor Deposition System, the CO2 Laser Sys- 
tem, and the analysis equipment. The Vapor Deposition System includes all of the 
equipment used to prepare the target: sander, polishing wheels, ultrasonic cleaner, 
pure-air system, vapor degreaser, deposit thickness monitor, and the VEECO VE-401 
automatic high vacuum evaporator. The CO2 laser system includes all of the equipment 
necessary to impinge a laser pulse onto the target: CO2 laser, target chamber, energy 
meter, and the optics. The analysis equipment includes all of the equipment used to an- 
alvze the results: optical microscope, and the scanning electron microscope. 
1. Vapor Deposition System 
The vapor deposition svstem consists of the equipment used to prepare the tar- 
get. [he equipments used to polish the target are the sander and polishing wheels. The 
equipment used to clean the target were the ultrasonic cleaner, the pure-air system. and 
the vapor degreaser. The film deposition was conducted using the VEECO VE-401 au- 
tomatic high vacuum evaporator and the deposit thickness monitor. 
a. Sander and Polishing Wheels 
To produce a smooth surface for film deposition the target underwent a 
standard polishing procedure using 240, 320, 400 and 600 grit sandpaper. This was fol- 
lowed by further polishing using the standard Buehler polisher with 1.0 micron and 0.05 
micron Alumina (a powder suspended in an aqueous solution). 
b. Ultrasonic Cleaning System 
The Somecor ultrasonic system (model SC-101TH) shown in Figure 9 was 
used to clean all of the targets. 
c. Pure-Air System 
In order to clean all particulates and dust off of the surface of the targets a 


clean environment is necessary. This clean environment was produced by the Model 
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Figure 9. Ultrasonic System 
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TM-48 pure air svstem made by the-Pure Air Corporation. This system filters out all 
QO.5 micron and larger particles and then blows the clean air into a small work area (00 
ein deep A 57 einen voll yocm wide): 
d. Vapor Degreaser 
This device consisted of a 1000 milliliter beaker, a target holder, a petri dish 
and a hot plate. A nuxture of acetone and ethyl alcohol 1s added and allowed to boil, 
the distilled vapors then condense on the target surface and run off by gravity. 
e. VEECO VE-401 Automatic High Vacuum Evaporator 
The VEECO VE-401 automatic high vacuum evaporator [Ref. 8] 1s shown 
in Figure 10. It consists of a vacuum system, pressure gauge controls and evaporator 


equipment. The vacuum system 1s composed of 


¢ mechanical pump capable of reducing the system pressure to approximately 10 
microns; 


e diffusion pump capable of lowering the systein pressure to the vicinity of 10-’ torr; 


¢ cold trap utilizing liquid nitrogen to remove condensible gases and prevent oil va- 
por from the diffusion pump from entering the bell jar. 
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Figure 10. VE-401 Automatic High Vacuum Evaporator. [Ref. 8] 
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The pressure gauge controls consist of 


¢ two standard thermocouple gauges capable of reading pressures from: 1000 microns 
of mercury tothe region of | micronvonincrcum 


¢ one thermocouple indicator-controller used to control the automatic functioning 
of the vacuum system, 


e one Bavard-Alpert tvpe ionization gauge capable of reading pressures from 
1 x 10°? tonpo 2 xe con, 
The evaporator equipment consists of various transformers and connections necessary 
to provide high current (up to 500 amperes) to the boat or filament used to hold and 
melt the substance to be deposited. 
ft. Deposit Thickness Monitor 

The SLOAN DITM-2A deposit thickness monitor [Ref 9) was usedmnme 
monitor the deposition process. The DIM-2A measures the mass of the deposit. The 
niass 1s measured by evaluating the change in frequency of a resonating quartz crystal 
as a deposit builds up on its surface. 

The DIM-2A consists of a sensing unit, an operating unit and a connecting 
cable. The sensing unit is a transistorized, crvstal controlled oscillator housed in a vac- 
uum tight package. The operating unit 1s an all solid state device whose principle display 
of information is in the form of a precision meter movement but also utilizes an audio 
amplifier and speaker principally for calibration purposes. | 

Before deposition the variable frequency oscillator is tuned so that the beat 
frequency between it and the crystal oscillator is zero (it is here that the audio amplifier 
and speaker are used in order to achieve more accuracy). After zero tuning the variable 
frequency oscillator frequency is kept fixed. As soon as deposition begins the crystal 
frequency changes and a difference frequency is heard. As the deposit thickness becomes 
appreciable (30 Angstrom units of Aluminum) the meter will begin to indicate. As the 
deposit thickness becomes greater the beat frequency will also rise according to the 


equauon 


f=AS> {32} 


Where 7 = deposit thickness in Angstroms 
AF = beat frequency in hertz 
p= bulk density of the depositing material in g/cy? 
A =a constant dependent upon the depositing material. 
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2. CO2 Laser System 
In order to produce a high intensity laser pulse, shape it and direct it onto a 
target under vacuum the following equipment was utilized: Lumonics CO2 laser, various 
optics, a target chamber, and an energy meter with its associated probe. 
a. CO2 Laser 
Weston [Ref. 10] described the CO2 laser in Figure 11 as follows: 


“A Lumonics TE 822 HP high energy TEA pulsed laser was used as the 
laser source. It is capable of producing 20 joules of output and has an adjustable 
pulse width from 0.05 usec to 5.0 usec.”” This laser uses an optically active medium 
of carbon dioxide, nitrogen, and helium, which is excited by a transverse electrical 
discharge. A population inversion results from the discharge electrons colliding with 
the CO2 and N2 molecules. Since one of the N2 molecular upper vibrational levels 
possesses an excitation energy which is approximately equal to that of the CO2 
molecule, the N2 serves to enhance the magnitude of the population inversion by 
the efficient transfer of vibrational energy during collisions between the CO2 and 
N2 molecules. Thus, the N2 serves as reservoir of energy which directly influences 
laser peak power and pulse length. 


“Helium is used to obtain a well diffused electrical discharge. The He 
also removes heat (generated by the electrical discharge) and lowers the gas mixture 
temperature. This results in less depletion of the CO2 upper vibrational level, and 
the interaction of the CO2 and He increases relaxation of the CO2 lower vibrational 
level. Therefore, the use of He allows the laser to be pumped at higher discharge 
currents, Which produces a larger population inversion and thus higher power. The 
high voltage power supply and gases were cooled by a Flowrite (RPCX mod 387) 
chill water unit with a flow rate of 25 gallons per hour at ..."[Ref. 10] 15 + 5°C. 


“The Lumonics laser with a multimode optic emits a 30 x 33 mm beam 
of approximately rectangular shape. The pulse shape can be varied in time from 
0.05 usec to 5.0 usec by adjusting the gas mixtures. For a variety of reasons, a 
standard long pulse mixture, resulting in a 4.55 usec pulse was used ....” ” The 
standard gas settings for long pulse are: 10 psig on the pressure gauges for CQO2, 
N2, and He, 8 standard cubic feet per hour (SCFH) for CO2 and N2, and 6 SCFH 
for He. This provided a reasonable pulse shape with a good distribution of the laser 
energy in time. The power density was varied by adjusting the voltage of the 
capacitor discharge ....”[Ref. 10] 


b. Target Chamber 
The target chamber is shown in Figure 12. As described by Wojtowich 
ieef. 11] it is 


“The VEECO 400 vacuum system... is utilized in conjunction with the 
CO2 laser for research of plasma surface interactions at the Naval Postgraduate 
School. The system is a modified vacuum deposition system which has been refitted 
with a vacuum chamber. The vacuum pumping system consists of a mechanical 
pump, a water cooled diffusion pump, and a liquid N2 cooled cold trap. The pres- 
sure range of the chamber is from 1.0 atmosphere down to 10-° atmospheres (atm). 
Three different gauges are required to determine the pressure in different ranges. 


74) 





Figure t1. Lumonics TE-822HP CQ2 Laser 


Pressures from one atmosphere (760 torr) down to 0.005 atmosphere (3.8 torr) are 
measured from a Matheson pressure gauge (model 63-5601) mounted on the top of 
the chamber; this gauge contains a Bourdon tube and socket.... Pressures between 
10-3 atmosphere (0.76 torr) and 10-‘atmosphere (7.6 x 10 “torr) are measuredsiaas 
thermocouple gauge located below the chamber. Pressures between 10-® atmosphere 
(7.6 x 10-6 torr) and 10-? atmosphere (5.5 x 10~ torr) are measured by an ionizamiem 
gauge located above the diffuston pump.” 


There were four ports mto the target chamber, two were of metal and one was a NaCl 
window. The laser beam entered through the fourth port which was of ZnSe. 
c. Energy Meter and Probe 
The energy meter and probe used were clearly described by Henson [Ref. 
12). 


“The energy output of the laser was measured using the Laser Precision 
Corporation Model Rj-7170 energy meter with the Model Rjp-736 energy probe 
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Figure 12. Target Chamber. 


which is designed to detect laser pulses of wavelengths between 0.35 and 11.0 
microns. The meter is capable of measuring the energy of a laser pulse of 1 nsec to 
I msec duration, and total energy ranging from 10 uJ to 10 J..... The energy detected 
at the probe is only a small percentage (0. 13°%) of the total energy of the laser pulse, 
whose flux on the target is calculated from the reflectance and transmittance of the 
ZnSe beam splitter, as well as the transmittances of the ZnSe focusing lens and ...” 
ZnSe window. 


d. Optics 
The optics used to direct and focus the laser light are as follows: 


® one 7.62 cm diameter ZnSe beam splitter with 99.38% reflectance, 0.13% 
transmittance, and 0.49% absorption at an incidence angle of 45 degrees; 


¢ one 7.62 cm diameter ZnSe focusing lens with 98.5% transmittance, 0.26% ab- 
sorption, and a 38 cm focal length; 


® one 7.62 cm diameter ZnSe window with 98.0% transmittance. 


De 


3. Analysis Equipment 

In order to analyze the surface of the targets the following equipments were 
uulized: optical microscope and a scanning electron microscope (SEM). 

a. Optical Microscope 

Preliminary analvsis of target damage and investigation of the concurrent 

pitting was conducted using a Zeiss ICM 405 Inverted Camera Microscope (Figure 13) 
with reflecting light. This microscope allowed viewing and photography of surfaces in 
magnification powers of 50X, LOOX, 200X and 500X. 





Figure 13. ZEISS ICM 405 Inverted Camera Microscope 


b. Scanning Electron Microscope 
Final target analysis was conducted using the Stereoscan 500 Scanning 
Electron Microscope [Ref. 13] which is capable of magnifications up to 280,000 X. Its 
basic function is to construct an image which has a three-dimensional appearance from 


the action of an electron beam scanning the surface. Under suitable conditions the re- 


solution of this instrument can be better than 4 nanometers with a depth of focus at least 
300 times greater than that of a light microscope at the limit of resolution. 
A basic diagram of the SEM is shown in Figure 14. The SEM is con- 
structed of five main components, these are 
e “Electron Gun produces a large, high intensity electron beam which is fitted into 


¢ the Column which controls and shapes the beam into a size useable for scanning 
microscopy. 


e The Scanning System scans the beam over the sample in a television type raster. 
The beam scanning over the sample releases electrons from it. 


e The Electron Collector Assembly collects these electrons and converts them to an 
image Which can be viewed by the operator. 


e Tlre Control Electronics allows vou to control the performance of items 1 to 
4."(Ref. 13] 
The onlv other capability of the SEM which was used for this thesis was its 
Elemental Analvsis function. 


“All samples will emit X-rays when struck by the electron beam. These 
X-ravs are characteristic of the element from which they are emitted. Bv the use of 
special detectors sensitive to X-ravs the sample may be analyzed for its constituent 
elements, Ihe results nay be either a plot of the concentration of the element 
present, or the spatial distribution of a chosen element on the sample surface. Par- 
ticles as small as 1 micron may be analvzed.”[Ref. 13] 


B. THE EXPERIMENTAL PROCEDURE 

This section explains the methods used to calibrate the equipment, prepare the tar- 
gets and conduct the shot. 

I. Calibration 

The only equipment item needing calibration by this researcher was the Deposit 
Thickness Monitor (DTM-3). Its calibration for copper is described below. 

Since no multiple beam interferometer was available, calibration for copper was 
achieved by measuring the initial mass (0.3130+.0001 g) and surface area 
(9.61 + .1cm7) of the crystal, then depositing to a beat frequency of 187.5+.25 kHz. The 
crvstal mass was then again measured and found to be 0.3265 + .0001 g. The mass dif- 


ference being .0135 g vields a film thickness of 


T = —mass 
area X p {33} 
a 


= 4.1 7S 
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Figure 14. Basic Schematic Diagram of the Operation of the SEM. [Ref 13] 
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Measurements on photographs of the crystal’s edge (see Figure 15) reveal a thickness 
of 


T=4.0um {34} 
Taking the average value of these two results vields a thickness of 
T = 4.089 p01 {o0} 
Manipulation of equation (32) thus yields the constant for copper to be 
mex i 


Be Ae {36} 
1.952 + .003 


A 


and so the equation for determination of the thickness of the copper deposits 1s 


| AF( Hz 
Loo x vaileae 
p(g/cnt ) 


with the density of copper being $.9Sg/cr. This equation 1s used throughout this thesis 


as the sole source for the determinauon of the film thicknesses. 
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2. Target Preparation 

The following procedure was developed using standard microscope slides. Films 
Were deposited onto these slides using various preparation and deposition techniques 
and then the film was analyzed. This analysis consisted of scrutiny under an optical 
microscope and also the placement of a light source on one side while looking for light 
leaks (holes) on the other. Correlation between these two methods of analysis yielded 
the general characteristics of film defects to be easily discernable from pitting caused by 
Unipolar Arcing. The techniques mentioned below yielded the best results, fewer than 
10 significant defects over an approximate area of 5 cnr, 

The target substrate was manufactured from one-half inch diameter type 304 
stainless steel rod which was cut into one centimeter sections. One face of these cylin- 
ders was polished to the extent obtainable with .05 micron Alumina using standard 
polishing techniques. 

In order to deposit a fine quality film it 1s necessary for the surface to Bervems 
cepuonally clean and free of all loose particles and dust (:nost of these particles have a 
much larger diameter than the thickness of the film). To obtain this necessary level of 
cleanliness the targets were first cleaned ultrasonically in an acetone bath. Further 
cleaning was done in a clean air environment produced by the Pure-Air system. In this 
clean environment the targets were dipped into acetone and then into ethvl alcohol, 
scrubbed with a standard drugstore variety toothbrush and dipped again. Since both the 
acetone and the ethyl alcohol contain particulate impurities the targets were also sus- 
pended above a boiling solution of acetone and ethyl alcohol in the Vapor Degreaser, 
thus the targets were cleaned further by distilled acetone and ethyl alcohol which con- 
densed onto their surfaces and ran off bv gravity. All transport of the targets after being 
in the Pure-Air system was conducted inside of a clean, sealed metallic box. 

After cleaning the targets were placed inside of the VEECO VE-401 Automatic 
High Vacuum Evaporator and allowed to outgas under high vacuum (approximately 
6.0 x 10-° mm Hg) for about one hour. The ‘irgets were placed 17 centimeters above 
the evaporation boat (hereafter simply referred to as boat). The DIM-3 oscillator was 
placed at the same height as the targets and equidistant from the point directly above 
the boat. The boat used was made of tungsten, it and its characteristics are shown in 
Figure 4, 

Before the 99.9% pure copper in the boat is heated a shield is placed over the 
boat to prevent deposition of impurities onto either the targets or the DI M-3 Oscillator. 


Current through the boat was slowly raised to 100 amperes. This raises the boat tem- 


perature to just above the melting point of the copper (1083 °C). After allowing all of 
the copper to become molten the current was further raised to 130 amperes. One minute 
more was allowed for the system to reach steadv state and allow any impurities or re- 
leased particles to settle. The shield was then removed allowing the vapors from the 
molten copper to deposit onto the targets and the DIM-3 Oscillator. The deposition 
was constantly monitored using the DIM-3 and when the desired change in frequency 
(and thus thickness) was achieved the shield was inserted and the current reduced to 
zero. A one hour cooldown under high vacuum (about 6.0 x 10-° mm IIg) was allowed 
before the target was removed. These films were then analyzed under the optical mi- 
croscope and their quality noted. 

The laser shot onto the target was conducted with the setup shown in Figure 
16. The target was placed under a vacuum in the 10-* mml{g range and shot with a fo- 
cused 4.5S5usec laser pulse. 

Final analvsis was conducted with the use of the SEM. The smallest pitting that 
penetrated through to the stainless steel and the largest pitting which did not penetrate 


was photographed as was their A-rav elemental analvses. 
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Figure 16. Equipment Configuration for the Laser Shot. 
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IV. EXPERIMENTAL RESULTS AND DISCUSSION 


This chapter reports the results of the experiments undertaken and gives an analysis 
of their implications. The beginning describes the results of experiments to determine 
the direction of emission of bulk material during the laser-target interaction. Structure 
was found in this emitted material by studving how it deposited onto a highly polished 
surface. This was followed by investigating damage to thin copper films on highly pol- 
ished stainless steel targets. Damage to these films led to a method of calculation of the 
arc current using basic thermodynamic principles. Finally, the results of experiments to 


find a relation between the size of the pits and their depth 1s given. 


A. INITIAL INVESTIGATIONS 

This section reports the results of experiments conducted to determine the direction 
of emission of bulk material during unipolar arcing. It also presents the results of de- 
position of this material onto a highly polished stainless steel target. 

I. Normal Emission of Bulk Material 

The first part of this research involved studv of the material ejected from the 
surface by unipolar arcing. Specifically, this section deals with finding out if there 1s any 
directional dependence to the emission of bulk material. 

Unless acted upon by external forces, such as electric and magnetic fields, a 
plasma, once produced, 1s expected to expand fairly uniformly to fill the vacuum and 
finally to deposit onto the surfaces of the chamber. To get a rough idea of the bulk 
distribution of the material ejected when the laser strikes the target the following exper- 
iment was conducted. A hole was made in a piece of aluminum foil and the foil placed 
normal to a copper target as shown in Figure 17. The hole was made to allow the laser 
pulse to penetrate and normally strike the target. The experiment was placed in a vac- 
uum of 2.0 x 10-* mmHg and shot with the laser 60 times at | Hz. Average pulse energy 
was 17.0 Joules. There was no visible copper deposited onto the foil. 

Investigating further it was found that keeping the foil parallel to the target 
surface but having the laser pulse come in at an angle (as in Figure 18) produced sig- 
nificant deposition. The parameters for this shot were the same as before, the only ex- 
ception being that the energy meter read an average pulse energy of 18.1 Joules. A 
qualitative examination of the deposition showed that most of the deposition occurred 


directly over the laser spot on the target. Different laser entrance angles produced the 
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Figure 17. Diagram of first experiment to determime the directionality of material 
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Figure 18. Diagram of second experiment to determine the directional dependence 


of material emission from a copper target: A= 1.7+0.2 cm. 
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same results. These experiments indicate that during the laser target interaction the bulk 
of the material which is emitted is emitted normal to the target surface. [t will be shown 
later that this should be universal for unipolar arcing. 

2. Deposition profile of ejected bulk material 

To further investigate this ejected material the aluminum foil was replaced with 
a highly polished type 304 stainless steel surface. The idea being that the ejected mate- 
rial, whether it was metallic vapors, plasma or molten liquid, would adhere to the sur- 
face. How it adhered would indicate the extent of the directionality and anv structure 
that may be present. 

In the first experiment the system was placed under a vacuum of 4 x 10-3 
mmHg. The laser target was then shot with 60 pulses at a frequency of 1 Hz. The av- 
erage energy of each pulse was 17.1 Joules, and the distance between the copper target 
and the steel surface, h, was 1.4+0.1 cm. Visual examination of the steel surface re- 
vealed the same results as with the aluminum foil. There was significant copper depo- 
sition normal to the laser spot of roughly the same size and shape (see Figure 19). 
Beyond this deposition region there was some spotted deposits of copper with the rest 
of the surface being covered with a blue film. This blue film indicates a deposition 
thickness of 100 to 200 nanometers over the rest of the surface. 

Examination under an optical microscope revealed several interesting phenom- 
ena. Figure 20 reveals significant pitting occurring in the region of dense copper depo- 
sition with the amount of pitting decreasing with the amount of deposition. This 
indicates that a hot and dense plasma is being emitted normally from the target surface 
and not expanding very rapidly in the radial direction. 

Figure 21 and Figure 22 reveal an intriguing ring structure found at two places 
on the edge of the region of dense deposition. This was also found in other shots to be 
described later. A possible explanation will be given then. 

Examination of the region outside of the dense deposition revealed circular 
copper deposits (Figure 23). There are no pits in these deposits. They are circular and 
appear to be fairly uniform. Each of these deposits are most likely caused by either the 
emission of a jet of copper vapor from a single unipolar arc occurring on the laser target 
surface or droplet formation (or some combination of the two). 

Reducing the vacuum to 5.3 +0.5 mmHg while leaving all of the other param- 
eters essentially constant produced very little deposition. However, it did produce sig- 
nificant pitting in the steel surface as shown in Figure 24, Figure 25and Figure 26. The 


phenomenon of rings formed around pits is also clearly evident. The figures show that 
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Figure 19. Steel surface showing bulk copper deposition. (6.4X\)( —— 1.0/0) 





Figure 20. Copper deposition on a steel surface showing a correlation between 
creater pitting density with greater copper deposition. (200N) 
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Figure 21. Ring structure: found surrounding a pit at the edge of the dense 





copper deposition region in Figure 19.(200X) ( Rj 50w1) 





Figure 22. Second ring structure: found surrounding a pit at the edge of the 


dense copper deposition region in Figure 19.(200X) (—@— 50us1) 






Figure 23. Circular copper deposits: found outside of the region of dense depo- 
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Fieure 24. Interference deposition pattern: found on a polished steel surface 


subjected to a copper plasma.(2U0N )) + 5071) 
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Figure 25. Second interference deposition pattern: found on a polished steel sur- 
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Figure 26. The region shown in Figure 25 at S50UX. ( ——4 20w71) 
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these rings form interference patterns when their central pits are close enough for them 
to interact. Figure 26 is a higher magnification view of Figure 25 and shows that these 
ring interference patterns are made up of many small dots of deposition. 

Further experiments did not vield such significant deposition as in Figure 19, 
however two deposition patterns were always evident: ring deposition and simple de- 
posits. Observations with the optical microscope reveal that every deposit 1s surrounded 
by a hght to dark blue film (as mentioned previously). This blue film is represented on 
these black-and-white photographs as the dark region surrounding the deposits. The 
data for Figures 27 through 33 are given in Table 1. In cach shot the target was hit 60 
times (at 1 E[z) by the CO2 laser. The results of shot C Were especially interesting ame 
surrounding blue regions of the deposits spread out in comet lke tails with all of their 


heads directed towards the same central region. 


Table 1. DATA FOR SHOTS A THROUGIEE: 
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These experiments show that the bulk of the material ejected from the laser spot 
is enutted normal to the target surface. While some of this material expands into the 
vacuum most of it comes out in jets. Pitting on the steel surface indicates that this ma- 
terial contains plasma but it also contains a significant amount of neutral metal vapors. 
The wave structure is most likely caused bv an oscillation in the incoming material which 
Was generated bv the central arc. Should this be the case, then the question of whether 
the rings are formed bv an outward propagating disturbance produced by a single arc 
or an oscillating arc arises. Also, if the velocity of these waves could be measured or 
calculated, then the period of the unipolar arc could be directly measured from the pe- 


riod of the rings. This would have to be taken up in later study. 
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Figure 27. Shot A: Centers are copper deposits surrounded bv yellow, and light 


and dark blue films on a copper background.(200N) ( 4 50) 





moure 23. Sitot Ay The center 1s a copper deposit. The copper rings dre on a 


blue background.(200N )( -—4_ SOun1) 
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Figure 29, Shot B: This shows copper rings on a dark blue background (dark 


roNon ea LOC — Se 





Fieure 30. Shot C: The center towards which the other deposits in thisssien 


point. (200N)( ——- S01) 
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Shot C: Ring deposition with the background pointing towards Fig- 


ure 30.(200N)( 4. 50m) 





Shot C: Ring deposition with its 


Figure 30.(50N)( H—4_ 0.27172) 
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Figure 33. Shot Drew These “Gre simple copper deposits with the usta) blue sie 
rounding film. the center is out of focus and steeply raised. (200X) 
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B. THIN FILMS AND THE UNIPOLAR ARC 

Placing a thin film of copper onto a lighly pohshed stainless steel target provides 
information on how the same unipolar are affects the two different metals. This section 
begins by providing data from these experiments further enforcing the argument that the 
primary surface damage mechanism 1s sunply local heating by the are current. Specil- 
ically, the unipolar are acts to melt and vaporize the surface imeisolaneee CClOr mee 
detailed theoretical description of the dynamucs of this process is given in Chapter Y. 
Lastly. this information, along with some basic thermodynamics, is used to estinate the 
arc current. The last section provides preliminary data for establishing aycorrelation 
between a pits depth and its width. 

1. Arc damage - a heating phenomenon 

To study the effect of the same laser shot. and the same un:polar arc, on two 

different materials a 0.37 micron copper film was placed over half of a type 304 stainless 
steel target. The target had been polished to the degree Obtainablewewan > 1icram 
Alumina. This target was shot ma vacuum of 3.5 x 10° mmllg. Optical and sean 


electron mucroscope analyses of damage in the laser spot revealed that while tlic copes 


film was penetrated, the copper appeared to have been pushed aside (Figure 34 to Fig- 
ure 37). This was even more apparent in the damaged areas outside of the laser spot. 
Figure 38 and Figure 39 show areas where the copper film was melted and somehow 
pulled radially outwards with no damage to the underlving steel. They also show regions 
where the underlving steel also suffered pitting vet the molten steel did not travel so far 
as did the molten copper. 

» These results are what would be expected given that the dynamics of the process 
is an electric arc joule heating and melting the material in a small region. This molten 
material then being acted upon by a radial force. The melting points of copper and iron 
feel) Grama tespectively. hus, for a given amount of heat placed over a 
short period of time, the copper would be expected to melt sooner and remain molten 
longer as it passed over a surface than would the iron (steel). 

2. Calculation of Unipolar Arc Current 
Assunung that the driving mechanism of the surface dvnamucs is produced by 
joule heating. then a rough calculation of the arc current is possible. This section will 
use some basic thermodynamics to develop a formula for this arc current. Since the 


currents of the unipolar arcs are most likely different as a function of time. this formula 
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Figure 34. Optical microscope image of a focused laser spot: 


with a 0.37 micron copper filim.(SOOX) (single shot)( FK_—- 20w) 
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SEM image of laser spot. 
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Figure 38. Optical microscope image of damage outside the laser spot: steel target 


overlaved with a 0.37 micron copper film. (SUON) ( -—Ae 20 purr) 
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Fieure 39, Optical microscope imag urther damage outside the laser 
Spot: steel target overlaved with a 0.37 mucron copper film. (300N) 
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is only meant to provide an order of magnitude estimate of the average are current. In 
light of this fact. several assumptions are made to simplifv the calculauons and the for- 
mula. Lustlyv. a calculation will be performed using an actual pit. 

Assuming that the pit is evlindrical with radius r and depth h. then the mass in 


the pit 1s 
: ) 
pis prroh {38} 


Where p is the density (8.95¢/cvr for copper). The amount of heat necessary to melt all 


of this and to vaporize one-quarter of it is 


Q, = (CAT, +L + CA T-0,252,). {39} 


+6 


Where AJ, = temperature change from initial temperature to melting point 
(1063°C for copper initiallv at 20°C) 
AT,;= temperature change from melting point to vaporization point 
C,= heat capacity for the solid (.092cal/g°C for copper) 
C,= heat capacity for the hquid (0.11 8cal/g°C for copper) 
L,= latent heat of fusion (49 cal g for copper ) 
L, = latent heat of vaporization (1250 cal/g for copper). 


Here the heat capacities are assumed to be constant over their respective regions. This 
assumption is fairly accurate over the regions given and should not be a significant 
source of error. Another assumption is that one-quarter of the material in the pit 1s 
vaporized. This is simply a very rough estimate which comes from observations as in 
Figure 39 and Figure 19. Figure 39 shows that a lot of material 1s going into the for- 
mation of the rim, whereas Figure 19 appears to show that a significant portion of the 
material in the pit is ejected. This assumption may be an order of magnitude source of 
error since the energy required to provide the latent of vaporization for all of the copper 
Dieetie pit is eight times that required to heat it to its boiling point 
(1118°C in a vacuum of 7.6 x 10-* mml 1g). The factor of one-quarter was chosen to 
account for both the observed bulk inaterial emission and the observed rim (which in- 
dicates that most of the material remained). The actual value of this factor is unknown. - 


In joule heating the power, P, provided by the arc current, I, 1s given by 
P=I°R. {41} 


R is the resistance over the current path and 1s given by 





R= Pr 4 {42} 


Pere p, 1S the resistivity (1./ x 10-*€2 «71 for copper). For simplicity the resistivity is 
assumed constant. 
Slee Porvcn ia tis Case 1S the rate Of production of enerey, and heat, Q,, is en- 


ergy, then 


Q, = et zs 
~ PR =) 
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Here t is the lifetime of the arc. It will be assumed that the current is constant 
throughout its lifetime and zero at any other time. A further assumption is that all of 
the heat produced by the arc goes into melting and vaporizing the metal in the pit. Of 
all of the assumptions in this derivation this is the one most subject to question. Cop- 
per’s high thermal conductivity assures that heat will be transmitted into the surrounding 
metal. Taking this heat loss into account would severely complicate this dérivation and 
so that must be left for further research. 


Finally, combining equations (42) and (43) provides the formula for the arc 





current 
ae 
<M 
j= /—-—. {4} 
hp, t 
The current density is then simply 
aa) 
| os ? 


{45} 





p,m ht 


These equations should be fairly accurate for anv arc in any material. 

Figure 40 displays the unipolar arc pit which will be used in the following cal- 
culation of an arc current. The film thickness, h. is 1.89 microns, the pit radmis, 1) isu 
microns. Given these parameters and making the assumption that the arc lifetime 1s one 


inicrosecond the heat required 1s 
O.=9 32 alla) {46} 
The arc current is then 


[=28.6 Amps. {47} 


This value is very near the current assumed by Schwirzke [Ref. 4, p. 611] and that re- 
quired by Wieckert [Ref. 14, p.500] providing further confirmation to their theories. The 


current density 1s 


ga Ole 10° Amps'cm’. {48} 
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M photograph of unipolar arc pit fully p 
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ing a 1.89 micron copper frlm. (S5000X) 


This value is an order of magnitude below that required by Werckert [Ref. 14. p.500]. 
This is most likely due to the copper having moved radially away from the central arc 
feeias beenpoOserved 1m the previous Section, 1-¢. a smaller cross section of the arc. This 
motion has the effect of increasing the value of the current density, however, since the 
amount of material involved in the calculation of the heat input is the same. the value 
eeeine arc Current remains the same. 

ie nietmourscttmainme the arc current will be much more accurate when some 
of the assumptions are removed. In particular, a good value for the are hfetime must 
be found. More detail in this calculation is meaningless without the hfetime. 

3. Pit depth experiment 

The purpose of this experiment was to find if there is a relationship between a 
pit’s depth and its diameter. To accomplish this copper films of thicknesses varving from 
One to ten microns Were placed onto polished tvpe 304 stainless stecl targets. These 
targets were then illununated with a single focused pulse of the CO2 laser under roughly 


the same conditrons. 
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Data was acquired using the scanning electron microscope. Its elemental anal- 
vsis function was used to scan for sharp differences in element concentration over the 
pit. An example of a pit which did penetrate fully through the film is given in 
Figure 41, Figure 42 and Figure 43. An example which demonstrates a pit which did 
not fully penetrate is given bv Figure 44 and Figure 45. The elemental analysis function 
will penetrate approximately one nucron into the surface. This is not considered to be 
a significant source of error since it is changes in the elemental concentrations which are 
looked for. The fact that all analyzed pits either showed a sharp change in concentration 
Or no apparent change at all further minizes this as a source of error. 

Figure 46 displays the results of this experiment. Raw data is provided in Ap- 
pendix A. The data show that there is no correlation between pit depth and diameter 
for pitting on the copper films. This conclusion comes from the fact that in many cases 
pits with relatively smaller diameter penetrated the film whereas other much larger ones 
did not. This conclusion may be specific to pitting on copper films due to the sharp 
differences between the shapes of these pits and those previously investigated on un- 
coated stainless steel. There is also dramatic differences between the pits on the same 
target. Figure 41, Figure 44, Figure 47 and Figure 48 are just a few examples. The 
large range 1n error in the power densities in Appendix A is due to the measurements 
of the small area of the laser spot. 

Figure 46 also shows that the maximum depth for pitting was less than 10 
microns. No pit was found which fully penetrated the 9.88 micron film. This includes 
pits inside the laser spot. Difficulty in acquiring time on the SEM prevented further in- 


vestigation. 
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Iron distribution of Figure 41. 


Figure 42. 





Figure 43. Copper distribution of Figure 41. 
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Figure 44. Example pit which does not penetrate the copper film. 
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Figure 45. Tron distribution in Figure TEs 
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Figure 46. Results of pit depth versus pit diameter experiment. 
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Figure 47. Example of Unipolar Arc pit diversity. (1320N) 


9,71KX 20KU WO'34mM 5: 68008 P: 98862 


 : [Eoeeennnceemeneemninenmenenemenencsenennnnamnnns 


Ory 
} A 
: a 


i 

UY VL ps 
Yay 

LP gi ph 

y% a f 

np i wy 





Figure 48. Example of Unipolar Arc pit diversity. (8710X) 
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V. THE SURFACE DYNAMICS OF UNIPOLAR ARCING 


The background theory combined with the experimental results provide an indi- 
cation of the dynamics involved in the surface interaction. Simply put, the surface sees 
a high concentration of ions producing an intense electric field over a sma!l region. This 
generates the emission of a lugh electron current which heats the surface producing 
melting and vaporization of material in that region. Indeed, this is quite similar to arc 
welding and so much of the physics of these two processes must be similar. This chapter 
will exploit that link. It begins with a theoretical description of the generation of the 
pit or kevhole as it is referred to in arc welding, and concludes with a theory on the 
formation of the rim. Much of the theory presented is specific to the macroscopic 
process of arc Welding and is only qualitatively related to the microscopic process of 


unipolar arcing. Further research will be necessary to make it specific to unipolar arcing. 


A. PIT FORMATION 

As mentioned in chapter II, the cathode spot is both emitting a high current of 
electrons and absorbing a large ion flux. Heat is mostly provided by the electron current 
with a small contribution from the energy of the 1ons (see Appendix B). The ion energy 
comes from their thermal energy, the energy gained in their acceleration through the 
cathode potential drop and their recombination energy. This recombination energy 1s 
approximately equal to the material’s ionization energy. All this energy serves to heat 
the small region of the cathode spot. [eating of the cathode spot gives rise to melting, 
vaporization and thermal electron emission. 

As described by Quigley [Ref. 15], in order for a pit to form the power density must 
be high enough for the combined effects of the ion pressure and the pressure caused by 
the vaporization of the material to produce a depression in the pool of molten liquid. 
Once formed the depth and shape of the pit will be determined by the forces involved, 
that is, the effects of those forces which tend to open the pit and those which tend to 
close it (see Figure 49). 

1. Forces which tend to open a pit 

The forces which tend to open and maintain the pit are the ion pressure (p,), 
vapor pressure (p,) and the recoil pressure (p,). In the laser-target interaction the radi- 
ation pressure would also have to be considered, but this will be neglected here since the 


plasma- surface interaction is all that is being described. 
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Figure 49. Forces involved in the formation of the pit. [Ref. 15, p.31]] 


The ion pressure comes froin the momentum imparted to surface by the impact 


of the ions. It is given by 


a 


Pp=> av (49} 
Where p, is the ion density and v is their velocity. Schwirzke [Ref. 4] gives a plasma 
density of 10'%c1-3? when one monolayer of atoms has been released from the surface. 
Schwirzke also stated an electron temperature of 18 eV. To find the upper limit on the 
beam pressure it is assumed that all of the ions enter into the sheath one debye length 
from the surface with negligible energy. The ions will be accelerated through the entire 


sheath potential. The velocity is then given by 


2eV 
is S 
v=,/ 7; {50} 


Where #7, is the ion mass. With p; = 71, the beam pressure becomes 





ii 








ie + pn 
HEN: 
7 2n; My {51} 
= eV 


1.45 x 10° Nm’. 


This is the upper bound on the beam pressure assuming a plasima density provided by 
the partial 1onization of one monolaver of atoms. Currently the value of the plasma 
density above the pit during arcing 1s unknown. 

The vapor pressure comes from the material vaporized inside of the pit. This 
pressure 1s dependent upon the temperature. Calculations of the temperature and thus 
the vapor pressure distribution inside of the pit are comphcated due to the dynamic na- 
ture of that environment. The vapor pressures as a function of temperature for alumi- 
num, copper and iron are given in Figure 50. Note the rapid decrease with temperature. 

Lastly, the recoil pressure 1s due to the vapors leaving the surface with a finite 
velocity and hence imparting a recoil pressure onto the liquid surface. According to 


Quigley {Ref. 15, p.317} the recoil pressure is given by 


Ww? 
p, = 


= ——— {5m 
A“p,Q° 


where W,A is the power density, p, is the vapor density and Q is the energy required to 
vaporize one kilogram of the material. From the calculation of the arc current in 
chapter IV W’A is 3.3 x 10 W/»r and Q is 1.94 x 10°) ‘ke. The vapor density canes 
found by taking the mass of the vapor ejected and dividing by its volume. Assuming that 
the surface temperature is 2200 kK, then the velocity of the vapors is 930 m's. If the arc 
lasts for one microsecond the vapors will travel 9.3 x 10m. Taking data from the pit 
used for the arc current calculation in chapter IV finds a vapor density of 4..*%kg'm’. 


The recoil pressure is then 
p, = 6.44 x 107 Non’, {53} 


The values presented here are only rough estimations accurate to only an order 
of magnitude. Nonetheless, it is clear that the recoil pressure is a significant force 


tending to create the pit in the hquified surface. The contribution of the ions depends 
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on their density. The ion density would have to be two orders of magnitude higher for 
the 10n pressure to become comparable with the recoil pressure. 
2. Forces which tend to close a pit 

The forces which will act to close a pit are gravity and surface tension. lhe 
gravitational pressure, pgh, at any point is simply equal to the weight of the metal above 
it. In Figure 49 his the depth at any particular point in the pit, and d 1s the maximum 
pit depth. This force should not be significant in pits formed by unipolar arcing. 

Surface tension can be thought of as the force which acts to minimize the sur- 
face areca of a liquid. In the case of pitting the surface tension force acts against the ion, 
vapor and recoil pressure forces as shown in Figure 49. Assuming that the surface at 


the bottom of the pit is a hemisphere, then the surface tension pressure is 


Be) 


Py a {54} 


where y is the surface tension coefficient and r is the radius of the hemisphere. The 
surface tension coefficient for copper varies from 1.03 to 1.15 N/m depending on the 
temperature (see Figure 51). Assuming a value of 1.1 Nm and a pit radius of one 


nucron gives a surface tension pressure of 


oy Se Nae {55} 


While this value is much greater that either the 10n pressure or the vapor pressure, it 1S 
much less than the recoil pressure 
3. Overview of the dynamics of pit formation 

The process of pit formation begins with the concentration of an electric field 
around a ‘whisker.’ This causes increased 10n bombardment of the whisker, heating it 
and increasing the emission of adsorbed gases and other neutral particles. Collisions 
with electrons ionize a small percentage (approximately 2%) of these neutrals which 
causes a dramatic rise over the background plasma density (about 4000X) above the 
Whisker [Ref. 4]. This rise in plasma density rapidly reduces the debye length with a 
subsequent increase in the electric field above the surface. When the electric field be- 
comes high enough arcing occurs. This is followed by a rapid temperature rise.in the 
cathode spot accompanied by melting and vaporization of the surface. This 
vaporization creates a recoil pressure which overcomes the surface tension pressure and 
pushes a pit into the surface. As long as the plasma density above the cathode spot re- 
mains high the pit will deepen. Since it 1s the high electron temperature along with the 
high local neutral concentration which generates the high local plasma density, the 
plasma density in the pit will decrease with pit depth. This is due to the plasma in the 
pit becoming isolated from the bulk plasma. When the electron concentration and en- 
ergv in the pit become sufficiently depleted such that ionization of neutrals is signif- 
icantly reduced, then the plasma density will drop. The subsequent reduction in the 
debye length lowers the electric field past the point where arcing stops. 

Andrews and Atthev [Ref. 17] give an excellent mathematical description of the 
generation of a pit, or ‘keyhole,’ by a continuous high power laser beam. Their pit is, 
however, on the order of a millimeter in diameter and the heat input 1s from the laser 
as it drills the hole. Further research may be able to modifv their solution to fit the 


physics of unipolar arcing. Their solutions to the hole profiles for different power den- 
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Figure 51. Surface tension coefficient versus temperature for copper. [Ref. 16] 


sities is shown in Figure 52. In this graph ‘a’ represents the beam radius. Plots A 


through E represent increasing power density. 


B. RIM FORMATION ; 

Two forces are responsible for the formation of the rim: the ejection of the vapor 
from the pit and surface tension. The vapor ejection serves to remove the liquid on the 
sides of the pit and pushes it radially outwards as it exits from the pit. The motion of 
the molten rim away from the pit will be shown to most likely be due to a gradient in 
the surface tension. 

I. Dynamics due to the ejected vapor 
The vapor produced in the pit exerts forces on the walls of the pit as shown in 
Figure 49. The rushing vapor produces a shear stress along the wall of the pit. Since 
that wall is a liquid-vapor interface the shear stress and the pressure gradient produce 
vertical flow which brings the liquid out of the pit.! At the exit of the pit the radial va- 


por pressure gradient pushes the liquid radially outwards to form a small rim around the 


pit. 


1 The liquid is also removed from the pit by the action of the recoil pressure displacing molten 
metal in the pit center. 
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Figure 52. Example pit profiles: equilibrium hole profiles calculated to be 
produced by a continuous high power laser beam whose diameter 1S on 


the order of a millimeter.[Ref. 17] 


The vertical flow also transports heat. This heat transport acts to produce more 
melting on the surface which explains the downward slope found between the rim and 
the pit (Figure 53 and Figure 54). 
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2. Dynamics due to surface tension 

Surface tension also acts to transport the surface liquid radially away from the 
region of the pit. How it does this depends upon whether or not the pit fully penetrated 
the material, in this case the copper film. If the pit fully penetrates it will be shown that 
the surface tension acts to reduce the surface area. In so doing the rim expands radially 
until the heat contained in the liquid is no longer sufficient to liquify more of the copper 
film. Two arguments will be given for the case where the arc does not fully penetrate 
the material of the film. The argument by Lancaster [Ref. 15, p.113] shows that when 
the surface tension decreases with temperature it will induce an outward radial flow on 
the surface of the liquid. Unfortunately, tvpographical errors in the reference produced 
a break in his argument. Attempts to resolve this produce an argument from which 
outward radial flow is produced by a surface tension which increases with temperature. 

Both arguments will be given concurrently. 

a. Full film penetration 

Since the pit depth is generally on the order of a few microns, then a film 
of that thickness or less on a surface will be easily penetrated by a unipolar arc. When 
this occurs the surface tension of the liquid acts to reduce its surface area. It will have 
a semicircular cross section and have a contact angle, 8, with the underlying surface 

(Figure 55). This contact amale sive py, 
cos0= Ree ' (oF 


VY 


1A 


where y5. 7, and y,, are the surface tensions of the underlying solid surface. the liquid 


film and the surface-liquid interface respectively (Figure 56). 
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Figure 55. Pit fully penetrating a film over a surface. 


The surface tension produces an outward radial force on the film 
surface(Figure 55). As heat is transferred from the liquid to the solid film more of the 
solid film will undergo the transition-to liquid and the ‘rin’ will progress radially away 
from its starting point. This will continue until there is no longer sufficient heat present 
for the transition to continue. At this point the heat will continue to flow out of the 
liquid and it will solidify in place while maintaining its shape. 

While this explains the damage found in Figure 39 it does not completcly 
explain that found in Figure 38. The large size of this damage indicates that it could 
not have been caused by a single arc in a small area. Thus either several arcs produced 
this damage or the arc followed the rim as it receded. The mechanism for the arc fol- 
lowing the receding rim may come from the reduced vaporization temperature for copper 
ina vacuum. The melting point for copper is 1083 °C. Dependmg on the vapor pres- 
sure, its boiling point can be very near its melting point. Ata pressure of 7.6 x 10-4 
mmHg the boiling point is 1118 °C. The melting point for the iron in the stainless stcel 
is 1535 °C. Thus the inner corner of the copper rim could have become the source of 
the vapors producing the plasma density needed for arcing. Thus the arc following the 
rim as it recedes may have been the heat source which generated the damage. 

b. Flow induced by a surface tension gradient 
The following argument deals with steady flow conditions in a Weld pool. 


While the process of unipolar arcing is on a inuch sinaller scale and is inherently dy- 
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Figure 56. Surface tension forces on a liquid drop in contact with a surface. [Ref. 
15, p.35] 


namic much of the same physics should be occurring. Thus this argument is not meant 
to present a full description ofa mechanism for liquid flow during umipolar arcing, rather 
it 1s only to present a ‘possible’ mechanism for this flow. More research is necessary to 
determine the validity of this arguinent in relation to the unipolar arc. Most of the fol- 
lowing argument comes directly from, Ref. 15. 

The general solution for-flow generated in a conducting fluid by an electric 
current Which diverges from a point source on the surface (Figure 57) was found by 
Lundquist (Ref. 19]. Lancaster (Ref. 15, p.113] apphed this to a pomt source at the or- 
igin entering a semi-infinite hquid whose surface is on the plane 6=7/2. The stream 


function, in spherical coordinates, for this case is given by 


ie 
Wj) = = [a)?+ bi +e+(142) In(1+2)] {57} 


ny 





whe.e uw, = magnetic permeability 
I = electric current 


yn = viscosity 
= coved 


and a, b andc are constants. The boundary conditions for this problem require that 


(1,0) = 0 eu). 
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i) — 0. {60} 


Figure 57. The current I enters at the origin of the coordinate system and flows 


radially into a conducting fluid. [Ref. 19, p. 90] 


When the svstem Is in steady state the net shear stress at the surface must 
be zero. Hence there is a balance between the stress due to the flow, p,, and that due 


to the surface tension. Thus 
e 3y 
ot or = =u (61) 


Lancaster then solves for the shear stress with the result that 





2 

| ”) | 

a ayy Teh dey) ingles i) | 2a+ ze i {62} 
nr? 1) 


Attempts to reproduce this result have failed. Calculations shown in Appendix C show 


the shear stress to be given by 


7 
fal 


pan te | Re [astsbieci4y In(1+A)]4+./ 1-27? (20+ £5 )} {63} 
yf 1-2 





An‘y? 
The difference lies in the numerator term under the square root sign in equation (61). 


Equation (62) shows that only the [—/? term in the numerator is under the square root. 


This is also the only term multiplied with the 2a+1/14/, and the sign in front of the 
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,/1—/? is plus in equation (62) while in equation (61) it is minus. It is most likely that 
tvpographical errors are the cause of much of the problem since for 4 = 0 equation (61) 
vields a complex value. However, when Lancaster combines equation (60) with equation 


(61) he gets 





Lol Cy 
— Toe (2a+1)+ ram 0 {64} 
so that 
An*r? aE 
+) aie (65) 
2 U 
Sinularly, combining equation (60) with equation (62) gives 
An’? a 
o--+| a {66} 
2 ly 


Where the subscripts L and D will refer to the different values founc by Lancaster and 
Downs respectively. 


The first boundary condition requires that c=0, the second requires that 








b=—a-—2\n2. The respective stream functions then become 
Cy OV 
12 7s en Anr? -- 
Va ee ey ee ee) 2 
LAR 5 2 om 2 me {67} 
ae . uo! | 


+(1+4) In(1+/)] 


a0 ae y 
NZ ugl*r l sual er 4] out l a Cr HARE 
— <— —_— <= at ee ae A” —— SE EeReTD = Ne Z 
i 4277 2 pale : Ul ic 


+(1+/) In(1+)]. 
In spherical polar coordinates the radial velocity is given by 


ATS A\ys 
Ve oton ee {69} 
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so that the velocities at the surface are 


i? oy 
pol” (3 ie 
[Milino =— —$—( 5--2in2-— (70} 
. An’nr \ 4 Hoa | 
Lu ie Qn’? ita 
? = O 3 Co 
dn’nr \ ¢ Hol 


In the absence of a surface tension gradient both cases vield a radial velocity 


given by 


ture 


ul? 
[Velaeo = - 5— { 


: 
axn"yr 


—2 In 2) {72} 


and the flow induced by the current 1s radially inwards on the surface. The physical 
significance between these arguments comes when there is a gradient in the surface ten- 
sion over the surface. Lancaster finds that the flow 1s reversed, 1.e. flow is radially out- 


wards, if the gradient in surface tension is finite and positive such that 


3 
Ta Oy 2 In 2. {73} 


Calculations shown thus far reveal that flow will be reversed only if the gradient in sur- 


face tension 1s finite and negative such that 





————— < -0.1137. {74} 


Given a point source of heat at the source where ©, is the power of the source, C, the 


specific heat and p the density, then the temperature distribution [Ref. 15] is 


gta eas 
- dapC,or 


Where o 1s the thermal diffusivity. The temperature gradient is then 
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dT GF 








—— = ee 76 
r dnpC,ar? _ 
If V is the power in volts, then Q, = VI and with 
cy oy aT 
cr ) ol mde i: 
then 
ny ak 
In’ ry? mele nV at 
ie ony {78} 
pele 2pCyouel - 
Lancaster's requirement for outward radial flow then becomes 
ee 2 
2pCyarpt| ° ig 
Whereas calculations show the requirement to be 
ac? 
nV oP 
eal a {SO} 


The surface tension for most pure elements decreases with increasing tem- 
perature, however smal] amounts of surface active elements such as oxygen and sulphur 
significantly reduce the surface tension (Figure 58 and Figure 59). As a function of 
temperature, high enough concentrations can lead to an increasing surface tension with 
temperature. Also, the change in surface tension with temperature of copper Is positive 
at lower temper.tures and negative at higher temperatures as shown in Figure 51. Thus 
either solution may generate outward radial flow. 

Lancaster also uses figures from an article by Matsunawa [Ref. 20] 
(Figure 60 and Figure 61 to prove his case. However, if Matsunawa 1s following the 
sign convention for flow in terms of a change in the stream function (Figure 62 on page 
73) then Figure 60and Figure 61 prove just the opposite. 


c figref ref figref page=no ref ** 
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Pure iron 
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Figure 58. Effect on the surface tension of iron by various non- metallic elements. 
(Reis 1S p27) 


Pure copper 
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Figure 59. Effect on the surface tension of copper by oxygen and sulphur. [Rel. 
ie. 39) 
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Figure 60. The velocity distribution, streamlines and isotherms: Marangoni flow 


with a negative gradient of surface tension with temperature [Ref. 15, 


pila 
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Figure 61. The velocity distribution, streamlines and isotherms: Marangoni flow 


with a positive gradient of surface tension with temperature [Ref. 15, 


p.118}. 
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Figure 62. Sign convention for flow in terms of a change in the stream function. 
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VI. CONCEGSIONS 


This research has provided a broad overview of the dynamics of unipolar arcing on 
the surface. It was shown that a significant amount of material is ejected during arcing, 
the bulk of which 1s ejected normal to the surface. This being independent of the inci- 
dence angle of the laser implies that it is the plasma produced by the laser pulse that 
caused the target damage. Ring structures found when the ejected material deposits on 
another surface indicates that the arcing oscillates and;or has a more complex structure 
than previously considered. Experiments also revealed no apparent connection between 
pit depth and diameter for the copper films. 

Experiments using thin films imply that damage increases as melting and boiling 
points decrease, however, this has not been rigorously shown. Other parameters may 
have significant contributions as well. The role of surface tension during ares which fully 
penetrate the film is clear. Bv reducing the surface area it produces a rim which expands 
radially. Radial expansion will continue so long as the molten metal contains sufhicient 
heat to liquifv the adjacent film. 

In invesugating the role of surface tension an apparent similarity between unipolar 
arcing and arc Welding was discovered. Exploiting this similarity has revealed a possible 
explaination of the dynamics of pit formation. The primarv force generating the pit 
appears to be the recoil pressure caused by the evaporating surface. This acts against 
the surface tensron of the molten liquid to push a pit in the metal. Molten metal 1s re- 
moved from the walls of the pit bv the action of the vapors as they rush out. The rim 
appears to be formed by the combined forces of the expanding gases (as they leave the 
pit) and radial motion caused by a gradient in the surface tension. It must be empha- 
sized that this theory is only preliminary. There are many differences between unipolar 
arcing and arc welding. These may mvalrdate much of this theory. 

Clearly, further research into the surface dynamics is required. The {mits ofitime 
research will provide better choices of materials for Tokamaks, spacecraft and any other 
devices which will come into contact with a plasma. Specifically, the following areas 


need to be further investigated: 


¢ how a material’s melting and boiling points affect damage and what role tts 
conductivity plays; 


e the cause(s) of the ring deposition patterns: 


¢ manv of the assumptions need to be removed and the theorv made specific to 
Unipolar Arcing. 
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APPENDIX A. RAW DATA FROM PIT DEPTH VERSUS DIAMETER 
EXPERIMENTS 


Table 2. DATA FOR PITS WHICH DID NOT 
PENETRATE THE COPPER FILMI. 


film thick- pit diam- vacuum laser pulse 
ness (um) eter (um) (10-* mmHg)! energy (J) 
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Table 3. DATA FOR PITS WHICH PENETRATED 
THE COPPER Tren. 


film thick- pit diam- laser pulse 
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APPENDIX B. ION ENERGY CONTRIBUTION 


The purpose of this appendix is to provide some good estimates of the energy con- 
tributed by the ions. The ions will contribute the energy they release as they combine 
with one of the surface electrons (hereafter referred to as recombination energy) and the 
kinetic energy thev acquire as they traverse the sheath. It will be shown that the re- 
combination energy 1s insignificant, whereas the kinetic energy is significant. For sum- 
plicity, throughout the following arguments it will be assumed that only the latent heat 
of vaporization 1s needed to vaporize the material. The ions must be able to provide this 


amount of energy if they are to reproduce their numbers. 


A. RECOMBINATION ENERGY CONTRIBUTION 

Kittel (Ref. 22. p.53] defines the cohesive energy, E,, as “the energy that musteee 
added to the crystal to separate its components into neutral free atoms at rest, at infinite 
separation, With the same electronic configuration.” For the purposes of the following 
argument the cohesive energy will be assumed to be the energy needed to vaporize the 


crystal. For copper this energy 1s 
E. = 3.49 eV, atom. {$1} 


The recombination energy provided by the ions will be approximately equal to the ma- 


terial’s ionization energy. For copper the recombination energy, E,, 1s then 
E,= 1-12 ex) lem {82} 


Thus the fraction of neutrals released that must be ionized so that the same number of 
neutrals can be vaporized bv the recombination energy alone 1s 
Es . 3:49 6N atom 


E, “ara ay ee a 0.45 ions/atom {$3} 





which is quite high. 

Since the majority of the ions (generated from the ejected neutrals) which return to 
the pit surface are produced in the sheath, then for simplicity it will be assumed that all 
of those ions are produced in the sheath. The number of neutrals which are ionized in 


the sheath is proportional to the ratio of the debve length, 4p, and the electron mean free 
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path, /,. The debye length was defined in chapter II, the electron mean free path 1s given 
by 


‘= {84} 


where n, is the neutral particle density and go, is their ionization cross section 
(~10-2m?). Given that each monolaver of atoms vaporized must undergo enough 
ionizations such that the ions produced deliver sufficient energy to the cathode surface 
to release another monolayer of atoms, then the neutral particle density can be taken to 
be 


n,=—— {$5} 


where p, is the particle density of one monolayer of the surface 
f=2 ilo x 10" jm? for copper). The fraction of neutrals which become iomized in the 


sheath is then 





Pp {86} 
=A/n i G; 


Thus recombination only provides .022/.45= .049=4.9% of the energy needed. This 
does not take into account the 1onizations produced by the electron current. This would 


entail a significant complication to this calculation and must be left for further study. 


B. KINETIC ENERGY CONTRIBUTION 

To find the maximum possible contribution of the kinetic energy it is assumed that 
the ions which go to the surface are produced one debve length from the surface. 
Schwirzke {Ref. 4] found that the vast majority of the ions generating the short debve 
length (and thus Unipolar Arcing) come from neutrals ejected from the surface. Thus 
it will be assumed that all of the ions are produced from neutrals evaporated from the 


surface. The kinetic energy that they provide upon striking the surface is then 


79 


my, {87} 


If the initial velocity is assumed to be negligible, then the ion velocity is given by 
v=al {88} 


where ais the acceleration and t 1s the trme for the ton to traverse the sheath. The ac- 


celeration ts given by 











eS 
lt Apm oa 
Given that the surface is at x=0, then the time, t, 15s 
ae (90} 
The maximum kinetic energy provided per ion is then 
KE =ely 
= ? ae n( ae ) — 


For an electron temperature of 18 eV, as given bv Schwirzke (Ref. 4 J, the kinetic energy 


of copper ions is 
KE = 88-4 eV, joi. {92} 


The fraction of the required cohesive energy which 1s provided by the maximum kinetic 


energy of the rons is then 


0.557. {93} 


NViN, is the ratio between the numbers of ions and neutrals as found tn the previous 
section (0.022). This estimates a significant contribution. 

[ons will be produced throughout the sheath, thus the actual contribution of their 
kinetic energy will be less than that given above. First it will be assumed that the ions 
are produced in a linear fashion such that one ion will be produced over each interval 


/,'N, Where .V, represents the total number of ions produced. It 1s also assumed that 
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the potential increases linearly with distance from the surface. The kinetic energy con- 


tribution from the ion produced in the jth interval from the surface 1s then 


: : J 
KE, = KE max “yo {94} 
The total kinetic energy contribution from all of the tons is then 
Ke . 
KE, = ) {95} 
j=} 
If f{x)=], then 
j+1 
| HICa Ges = {96} 
/ 
so that 
N; N; 
inet ay 
» | Weigee | fiwjdx = Dj. (97} 
j=l j=) 
The upper bound on the total kinetic energy 1s then 
omer gial 
KE,,, =——— | fix)dx 
Nj 1 
] I ve) 
oe N;KEmax— ON, KK Emax: 





ive 
=o K Emax: 2) 


Thus the ion’s kinetic energy provides approximately 27.9% of the total energy required 
for vaporization. 
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APPENDIX C. DETAILED CALCULATION OF THE SHEAKSii2-=] 


Whenever there is a viscous fluid moving with a nonuniform velocity shear stresses 
will be present. These stresses can be considered to act as if their is a frictional force 
between the elements of the fluid. This force coming from the ‘friction’ generated be- 
tween adjacent elements moving at different velocities. Figure A-1 shows the viscous 
stress components in the spherical coordinate system (1r,6,@). These components, as 


derived by Lancaster [Ref 15. p.Solrare 
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On a planar surface defined by @ = 2.2 the stress which acts to counter the surface ten- 
sion 18 p,. The purpose of this appendix is to calculate this stress component. 


Equation (56) gives the stream function to be 





ie 
W(r4) = fal?+bs-+e+(L- 2) In(1+4)] (56} 


4x77 
where 4 =cos@. The first boundry condition 
¥(r,0) = 0 {57} 


unplies that c=0. In order to calculate the shear stress the velocities in the radial and 


theta directions are required, these are provided by the stream function as follows: 


CO 
ty 





Figure 63. Components of viscous in the spherical coordinate system. [Ref. IS, 
p.85] 
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The required partial derivatives are 
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=—,/1-/ a {105} 











and 








a V [ 
i (2 a EE ae Pov ie yee {106} 


Placing these equations into the equation for the shear stress yields 
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